General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



NASA 

Technical 

Memorandum 


NASA TH -86515 


WE I BULL DISTRIBUTION BASED ON MAXIMUM LIKELIHOOD 
WITH INTERVAL INSPECTION DATA 


By Mario H. Rheinfurth 


Systems Dynamics Laboratory 
Science and Engineering Directorate t 


(*, SA-Tfl-86515) HI HULL DISINIEOTIOM EASID 
ON HAIIHUfl L1KE1IBOOD HIIH IN1EBVAL 
I NSFECTICS DAIA (IASA) 1 I p EC A02/HF A01 

CSCL 12A 

. 63/6 5 


July 1985 


M85-32647 


Uiiclas 

21999 


WASA 

National Aeronautics and 
Space Administration 

George C. Marshall Space Flight Center 



MS^C • Form 3190 <Hav. May 1983) 


TECHNIC Al 

. REPORT STANDARD TITLE PAGE 

1. HE PONT NO. 

NASA TM-86515 

2. GOVERNMENT ACCESSION NO. 

3. RECIPIENT'S CATALOG NO. 

4. title ano sustitle 

Weibull Distribution Based on Maximum Likelihood 

5. REPORT DATE 

Julv 30. 1985 

with Interval Inspection Data 

6. performing organization cooe 
ED01 

7. AUTHOR! S) 

Mario H. Rheinfurth 

8. PERFORMING ORGANIZATION REPORT tt 

9. PERFORMING ORGANIZATION NAME AND ADDRESS 

George C. Marshall Space Flight Center 

10. WORK UNIT NO. 

Marshall Space Flight Center 

, AL 35812 

11. CONTRACT OR GRANT NO. 



13. TYPE OF REPOR'i 4 PERIOO COVERED 

12. SPONSORING AGENCY NAME ANO AOORESS 

National Aeronautics and Space Administration 
Washington, DC 20546 

Technical Memorandum 



14. SPONSORING AGENCY CODE 

J 

15. SUPPLEMENTARY NOTES 

Prepared by Systems Dynamics Laboratory 

J 

It. ABSTRACT 1 


This technical note determines the two Wei bull parameters based upon the method 
of maximum likelihood as presented in memorandum, "Oxidizer Turbine (HPOTP) First 
Stage Blade ReMability Analysis," dated July 10, 1985, by the author. The test 
data used were failures observed at inspection intervals. The application was 
the reliability analysis of the SSME oxidizer turbine blades. 


i 


i 


17. KEY WORDS 

Reliability, Weibull Distribution, 
Maximum Likelihood, Interval Data 

18. DISTRIBUTION STATEMENT 

Unclassified-Unlimited 


19. SECURITY CLACSIF. <« t thla report! 

20. SECURITY CLASSIF. (of th'.e page) 

21. NO. OF PAGES 

22. PRICE 

UNCLASSIFIED 

UNCLASSIFIED 

10 

NTIS 


tare - For* SIM (Map 1M9) 

for tale by National Technical Information Service. Springfield. Virginia 221 61 


TABLE OF CONTENTS 


Title Page 

INTRODUCTION 1 

DISCUSSION 1 

NUMERICAL RESULTS . 4 

A. Unchamfered Blades (Table 1) 4 

B. Chamfered Blades (Table 2) 5 

CONCLUSIONS 5 

TABLE 1 . Unchamfered Blade Test Data 6 

TABLE 2. Chamfered Blade Test Data 7 

SIGNATURE PAGE 8 


iii 


PRECEDING PAGE BLANK NOT FILMED 


WEI BULL DISTRIBUTION BASED ON MAXIMUM LIKELIHOOD 
WITH INTERVAL INSPECTION DATA 


Mario H. Rheinfurth 
Systems Dynamics Laboratory 
George C. Marshall Space Flight Center 
Huntsville, Alabama 35812 


INTRODUCTION 


This technical note determines the two Weibull parameters 
based upon the method of maximum likelihood as presented in 
memorandum, "Oxidizer Turbine (HPOTP) First Stage Blade 
Reliability Analysis, " dated July 10, 1985, by the author. 
The test data used were failures observed at inspection 
intervals. The application was the reliability analysis of 
the SSME oxidizer turbine blades. 


DI SCUSSI ON 

This memo presents the results of a reliability analysis of 
the HPOTP first stage blades. The failure mode investigated 
here is a shank cracking of the blades on the downstream pres- 
sure side. It was assumed that the time distribution of this 
failure mode follows the two-parameter Weibull distribution 
whose probability density function is given by 

(1) f(t) - (B/n B )t 0-1 EXP[-(t/n) 0 ]. 


The parameter 8 is called the shape parameter and is positive. 
The parameter n is called scale parameter and is also posi- 
tive. Because the scale parameter n always represents the 
63.2 percent point, it is also called the "characteristic 
life". It has the same units as t, namely, time, whereas B is 
a dimensionless number. 


The Welbull reliability function (survivorship function) is 

( 2 ) R(t) > EXP - [(t/n) 8 ]. 

This is the probability of a blade surviving beyond time t 

Two sets of teat data were analyzed, one set for unchamie^ed 
blades (Table 1) and one set for chamfered/blended blades 
(Table 2). 

Estimation of the two Welbull parameters 0 and n was based 
upon the method of maximum likelihood. This method was 

selected for two reasons. First, it is capable of using all 
the statistical information contained in the test data. 
Second, the method utilizes the intuitively appealing princi- 
ple to determine the two parameters in such a way that the 
sequence of observations that actually occurred is the one 
having maximum probability. Usually the observations are a 
random sample r.f independent observations from the same 
distribution. Most reliability analyses assume that each 
failure time of a component is known exactly. However, for 
the present test data only the intervals in which failures 
occurred are known. The likelihood function can be set up to 
properly account for this condition. 

Suppose that R blades have been found to fail in the interval 

t t . < t. < t 2 . (i-1, 2 R) and S blades have been found to 

survive 1 beyond time t, (J - 1, 2, ...S). Then the sample 
likelihood function is given by 

(3) L = n[EXP-(t 1 ./ n ) 6 - EXP - (t 2 ./r,) 6 ] x l EXP - (t 3i /n) 8 . 


where tl.e first product n is the probability of R blades fail- 
ing within their respective test intervals and the second 
product ^ is the probability of S blades surviving their 
respective test runs. 

For the following analysis, it is convenient to work with the 
logarithm of the above likelihood function, which is 

R S v 

In L = l In [EXP-(t x ./n) e - EXP - (t 2 ./n)”] - (l/n) B f t 3j . 


( 4 ) 


OUJGINAL PAG*. ,S 
Op POOR QUALITY 


The given test data represent a special case in which the 
upper failure interval limit is equal to the total test time, 
i.e., t 2 . - t,.. Moreover, there are M-78 blades per disk and 
it is observedat the end of the i-th test that R. blades have 
failed 3ince the inspection time t t . and that t M - R . ) blades 
have survived the test. Supposing that there are a total of 
N test runs, then eq (4) can be modified accordingly to yield 

N N 
(5) In L = l R. MEXP-U^/n) 0 -EXP-(t 2 ./n) 8 ] - (1/n) 8 J (M-R.)t 8 . . 


Several numerical methods exist to maximize the log-likelihood 
function of eq (5). The results cf this analysis were 
obtained by the Newton-Raphson Method which finds the maximum 
of the log-likelihood function by solving the set of equations 


( 6 ) 


3 in L 

0 , 

3n 


3 In l 

= 0 . 

36 


This was done on a programmable hand calculator. The method 
works well if the starting parameters are close to the maximum 
likelihood estimators. With some experience, it is not very 
difficult to find good starting values. 

As mentioned earlier, the reliability function of eq (2) rep- 
resents the probability that a single blade will operate 
without a crack for at least t seconds. Of practical impor- 
tance is also the reliability of a turbine disk, which for the 
HPOTP first stage contains n*?8 blades. Considering the tur- 
bine disk to be a system composed of n identical and 
independent components the disk reliability is given by the 
cumulative binomial distribution as: 

(7) R„(k. t) * I (") R n " x (1 - R) x . 

where R is the reliability of the single blade. The reliabil- 
ity R D ( k , t ) is the probability of not more than k components 
falling within the time interval from 0 to t. This "system" 
reliability is, of course, always smaller than the reliability 
of a single blade. 


3 


ORIGINAL PAv& ^ 
OF POOR QUALITY 


Another statistical quantity of interest is the conditional 
reliability. This is the probability of a component (or sys- 
tem) to survive a specified operational time T after it nas 
already been lr op ration for some known time T . For the 
Weibull distribution, it is given by 

EXP - [(T + T 0 )/n] B 

(8) R r (T/T 0 ) = 

L EXP - 


(T 0 /u)* 


NUMERIC AL RESULTS 


Only the estimates of the Weibull parameters and the Bi-life 
will be presented. The Bi-life represents the operating time 
for which the reliability is 99$ (1)1 unreliability). The 
shape parameter 8 may be used as an indicator of the underly- 
ing failure mechanism. If 8 < 1 the failure rate decreases 
with the age of the component. This is often referred to as 
infant mortality or as burn-in failure mode. For 8*1, the 
Weibull distribution is identical to the well-known exponen- 
tial distribution, which represents a random failure mode. In 
this mode, the failure rate is constant which means the compo- 
nent does not age. This condition is more appropriate for 
complex systems whose components have different ages due to 
previous failures and replacements. A shape parameter 8>1 
indicates a wear-out failure mechanism. However, it is impor- 
tant to understand that a more accurate assessment of the 
underlying failure mechanism has to be based upon the confi- 
dence intervals of the Weibull parameters. The character i st 1 c 
life n is, in general, not of interest because it corresponds 
to a reliability of only 36.8$. When comparing different sets 
of Weibull parameters, it will be noticed that even though 
their values may differ substantially, the corresponding 
Bl-lifes or B.1 lifes, that is high range reliabilities, will 
differ only slightly. 

A. Unchamfered Blades (Table 1) 


Scale parameter 
Shape parameter 


B1 -life 


n - 1,142 ,538 sec . 
8 * 0.581 
t , - 4 1 3 sec . 
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B. Chamfered Blades (Table 2) 


1. Blended blades only (exclusion of data set 

marked by asterisk). 

Scale parameter n - 69,32*1 sec. 

Shape parameter 6 - 1.806 

Bi-life tj = 5.^29 sec. 

2. All chamfered blades (inclusion of data set 

marked by asterisk). 

Scale parameter n * 898,9*19 sec. 

Shape parameter 6 * 0.8208 

Bi-life t, - 3,309 sec. 

CONCLUSIONS 

1) The given test data indicate a significant statisti- 
cal difference between the chamfered and the unchamfered 
blades, the former showing an order of magnitude higher 
Bi-life than the former. 

2) There appears to be also a significant change when 

the single test run having unblended blades is included in the 

analysis. Especially worthy of notice here is the change of 

the shape paramotar 6 from a wear-out failure mechanism (B>1) 

to an infant mortality condition. This could be indicative of 
a quality control problem. Whether this change can be attrib- 
uted to the difference between blended and unblended blade 
design or to some other factor cannot be determined since only 
this single test run of unblended blades was available. 

With the above given Weibull parameters it requires only a 
small computational effort to determine single blade and disk 
reliabilities for various time points of interest using the 
appropriate equations presented earlier. 


TABLE 1. UNCHAMFERED BLADE TEST DATA 
(No Inspection time available.) 


Number of 

Cracked Blades Total Test Time (Sec ) 


2 

7 

1 

0 

1 

0 

5 

3 

13 

1 

0 

4 
2 

6 
0 
0 
0 
0 


10155 

6474 

6092 

3100 

4298 

3000 

3226 

3152 

3095 

3048 

2850 

2406 

2373 

2180 

1761 

1615 

1343 

721 
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TAME 2. CHAMFERED BLADE TEST DATA 


A. TCSTS WITH CRACKED BLADES 


No. of Cracked Biades 

3 

1 

1 

1 

1 

3 

1 

13* 


Inspection Time (Se c) 

4000 

2100 

1100 

1200 

1100 

0 

0 

0 * 


Total Test Time (Sec) 

5996 

4252 

3441 

3360 

3281 

3112 

1263 

1798* 


•Represents unblended blades. 


B. TESTS WITHOUT CRACKED BLADES (No inspection time necessary) 


Total Test Time 
(Sec) 

Total Test Time 
(Sec) 

7840 

2448 

6332 

1847 

5982 

1508.5 

4653 

1427 

4391 

1260 

4135 

1015 

3491 

895 

3241 

766 

3155 

756 

3077 

750 

2946 

666 

2880 

501 

2868 

371 

2849 

300 

2810 
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